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Abstract: We provide a hydro dynamical description of a holographic theory with bro¬ 
ken translation invariance. We use the fluid/gravity correspondence to systematically 
obtain both the constitutive relations for the currents and the Ward identity for mo¬ 
mentum relaxation in a derivative expansion. Beyond leading order in the strength of 
momentum relaxation, our results differ from a model previously proposed by Hartnoll 
et al. As an application of these techniques we consider charge and heat transport in 
the boundary theory. We derive the low frequency thermoelectric transport coefficients 
of the holographic theory from the linearised hydrodynamics. 
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1 Introduction 

Perhaps the most useful and widely studied application of holography has been to 
understanding the transport properties of strongly coupled systems [1, 2], In general, 
one is faced with a paucity of techniques to calculate real-time transport in theories 
that do not admit a quasiparticle description. In contrast, holography provides a simple 
prescription in which one can extract the linear response coefficients of certain strongly 
coupled held theories from the classical perturbation equations of black holes. 

The ultimate hope is that the ability to perform these calculations might yield new 
physical insights into transport at strong coupling. Particular targets for where such 
insight is necessary include the transport properties of the strange-metals. It is ex¬ 
tremely difficult to reconcile the experimental phenomenology of such materials with 
a quasiparticle or Drude picture of transport [3]. As such, it is important to develop 
techniques to study more general models of heat and charge transport. 

Motivated by this goal, there has recently been a large amount of progress in calculat¬ 
ing the transport properties of holographic models with broken translational invariance 
[4-13]. In particular, for a large class of theories it is now possible to obtain analytic 
expressions for all the DC conductivities in terms of horizon data [9, 10, 14-18]. More 
recently, the low frequency AC conductivity has been calculated in a simple example 
by exploiting a clever decoupling of the bulk currents [19]. 
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However, despite this progress, it remains unclear what the underlying physical pro¬ 
cesses that govern transport in these holographic models are. For translationally in¬ 
variant theories, a physical understanding of the boundary theory, at least at high 
temperatures, is provided by relativistic hydrodynamics. This allows a systematic ex¬ 
ploration of the physics of the boundary theory in terms of a derivative expansion of 
the hydrodynamic helds [20, 21]. 

The dehning data of relativistic hydrodynamics are a set of constitutive relations 
which express the energy-momentum tensor and electrical current in terms of 
a local fluid velocity, charge density and energy density. The dynamics of these helds 
can then be determined by solving the conservations laws for the currents 

= 0 

= 0 ( 1 . 1 ) 

The connection between the bulk geometry and the huid-dynamics of the boundary 
theory has been developed in a beautiful series of works (especially relevant are [22- 
25]). In particular, by performing a derivative expansion in the charge density p(x), 
energy density e{x) and huid velocity u^{x) it is possible to derive both the constitutive 
relations and the conservation equations of the boundary hydrodynamics. 

In this paper we use this connection to systematically derive a hydrodynamical de¬ 
scription of the simplest holographic model with momentum relaxation. Specihcally, 
we consider boundary theories where translational invariance is broken by turning on 
linear sources 0^^ = kxj\^ for massless scalar helds 0_4 in the bulk [11]. The breaking 
of translational invariance modihes the conservation equations of relativistic hydrody¬ 
namics to 


d,r = 0 

= ( 1 . 2 ) 

where the Ward identity for the stress tensor controls how momentum relaxes to equi¬ 
librium through scattering oh the scalars. By performing a derivative expansion both 
in the huid variables and the sources we can extract the boundary theory hydro¬ 
dynamics. This allows us to obtain the constitutive relations for J^, and {O up 
to second order in our expansion. 

Armed with these constitutive relations we can discuss the linear response of the 
boundary theory. At leading order in the derivative expansion our model agrees with 
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a previous proposal of Hartnoll et al [2] to include momentum relaxation simply by 
modifying the conservation law for momentum to 

where is a phenomenological momentum dissipation rate that we calculate. How¬ 
ever, beyond leading order, we will see that this proposal is too simple to describe the 
holographic models. The presence of the scalar held introduces corrections both to 
the constitutive relations and to (1.3) that are crucial to obtain the correct transport 
properties. 

By carefully taking these factors into account, we calculate the thermoelectric re¬ 
sponse coefficients of the boundary theory perturbatively in k. Our results reproduce 
the usual formulae expressing the DC transport coefficients in terms of horizon data. In 
addition, it is straightforward to extract the low-frequency AC response. This decom¬ 
poses into the sum of a coherent piece, arising from the part of the currents proportional 
to the huid velocity, and an incoherent contribution arising from the derivatives of the 
hydrodynamical helds^. In short we have constructed, for the hrst time, a hydrodynam- 
ical description that reproduces the known features of transport in these holographic 
models. 

The rest of this paper is organised as follows. In Section 2 we dehne our bulk model 
and explain briefly how to apply the fluid/gravity correspondence to discuss momentum 
relaxation. Our main results are the expressions for the constitutive relations of the 
energy momentum tensor (2.32), electrical current (2.35) and scalar expectation values 
(2.41). In Section 3 we linearise the hydrodynamics around equilibrium and calculate 
the thermoelectric transport coefficients of the boundary theory. Finally we end with a 
brief discussion of the possible future applications of these techniques in Section 4. An 
appendix is provided to describe the technical details of our fluid/gravity calculations. 

2 The fluid/gravity correspondence 

In this section we wish to describe how to apply the fluid/gravity correspondence 
in the context of the simplest holographic model in which we can study momentum 
relaxation [20]. We will therefore use the Einstein-Maxwell action in 5 bulk dimen¬ 
sions, supplemented by a set of massless scalar helds which will be used to break the 

^The same decomposition was recently found in a related model using completely different tech¬ 
niques in [19]. 
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translational invariance of the bonndary theory^ 


S 


1 

IGtiGm 



R + 12- - ^9'^'^dM>l>Adr,4>A 


( 2 . 1 ) 


Here capital Roman indices M, N refer to the bnlk coordinates whilst the index A = 
1,2,3 rnns over the spatial directions of the bonndary theory. The motivation for 
inclnding several scalar helds is that this allows the constrnction of bnlk solntions 
which break translational invariance in the bonndary bnt preserve homogeneity and 
isotropy in the bnlk. In particnlar there exist explicit black hole solntions to (2.1) 
when one applies the sonrces 0^^ = kxj, in the bonndary. Varions aspects of the 
thermodynamic and linear response coefficients of these solntions have been intensely 
stndied in [11, 19, 26]. 


The goal of this work is to provide a hydrodynamical description of the bonnd¬ 
ary physics dnal to (2.1). Onr starting point is the familiar translationally invariant 
Reissner-Nordstrom black brane. If we denote the bnlk coordinates by (r, x^), with r 
the holographic radial direction, we can write this black hole metric as 


gMNdx^ dx^ 
AMdx^ 
(pA 


—2u^dx^dr — f{r)u^Uydx^dx'^ -|- r'^P^^dx^dx'^ 


\/3g 

2A 


u^dx^ 


( 2 . 2 ) 


where we have introdnced a constant 4-velocity aligned in the ingoing-Eddington Finkel- 
stein direction, i.e. u^dx^ = —dv, and = rj^y + u^Uy projects onto directions 
perpendicnlar to u^. 

Notice that, becanse they are massless, we can tnrn on constant modes of the scalars, 
0 ( 4 , withont affecting the backgronnd geometry. As snch, /(r) is jnst the familiar 
emblackening factor of the Reissner-Nordstrom black brane 

f{r) = l — b/A + q'^/r^ (2.3) 


which has a horizon at radins tq satisfying 

( 2 . 


^Note that the fact we are working in 3-1-1 boundary dimensions is for ease of comparison with the 
majority of the seminal fluid-gravity literature [22-25]. We expect that the results presented here will 
generalise straightforwardly to, for instance, 2-1-1 dimensional theories. 


b = 4 


1 + 
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The thermodynamics of the boundary theory is well known. The chemical potential, /i, 
charge density, p, energy density e, pressure P, temperature, T, and entropy density, 
s, are given by 


^ ^ ^ q ^ ^ 3b p ^ h 2" = —1^2 - — 

^ 2rl ^ AttGn IQttGn IQttGn 2,7i \ r® 

Finally the expectation values of the scalars are trivial 


s = 


r 


3 

0 


AGn 


{Oa) = 0 (2.5) 

The derivative expansion 

We now wish to construct the hydrodynamics of the boundary theory in a derivative 
expansion. This correspondence, pioneered in [22], provides an algorithmic way to 
construct the hydrodynamics dual to a bulk action and has been subsequently applied 
to many different models. In particular, the fluid gravity correspondence has been 
developed in detail for the Einstein-Scalar action in the context of driven fluid flows 
[25]. Obtaining the hydrodynamics dual to (2.1) essentially amounts to generalising 
these calculations to a charged geometry [23, 24]. 

Whilst the application of the fluid/gravity correspondence is conceptually straight¬ 
forward it is rather computationally intense. In this paper, our main focus is on un¬ 
derstanding the hydrodynamics of the boundary theory. As such, we will be somewhat 
schematic in presenting the technical details of the gravity calculations. The specihc 
computations we have performed are explained in more detail in Appendix A. 


The starting point of this algorithm is to promote the parameters of the background 
solution u^, q, b, (j)^ to functions of the boundary coordinates to obtain 

= —2Ufj_{x°')dx^dr — r^/(r, q{x°‘), b{x^))u^{x°‘)uy{x^)dx^dx'' -|- r^P^y{x°')dx^dx'' 


/i(o)^ M \/3q'(x“) , , 

A\Jdx = — — — u^[x°')dx^ 

A = ‘Pa(x°) 


( 2 , 6 ) 


where we can think of M^(a;") as corresponding to a local fluid velocity in the boundary 
theory satisfying = —1. Now in general (2.6) is not a solution to the equations 
of motion of (2.1). Nevertheless, the essential idea is that, because it satishes the 
equations of motion on constant solutions, the corrections to (2.6) must be proportional 
to derivatives of the helds M^(a:“), g(a:“), 5((r"), 0^((r“). These correction pieces can be 
determined by solving the Einstein-Maxwell-Scalar equations order by order in this 
derivative expansion. 
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Before we begin solving for these correction pieces, there are two aspects of onr 
derivative expansion we shonld highlight. Firstly, we wish to implement momentnm 
relaxation in the boundary theory through position dependent sources for the scalar 
helds. Just as for constructing black hole solutions, the discussion is simplest if we 
choose the isotropic sources 

(j)^ = (2.7) 

These sources (2.7) dehne the laboratory frame of our theory - we will continue to use 
the calligraphic index A to refer to the spatial directions in this frame. 

Note, however, that the stress tensor of the massless scalars is quadratic in these 
derivatives. As such, these sources do not influence the metric and gauge held, and 
hence the boundary constitutive relations, until order (9^0^^)^ in the derivative ex¬ 
pansion. Physically this manifests itself in the fact that at leading order momentum 
will relax in our models at a rate ~ k"^. To see the effects of this scattering it is 
therefore necessary to study huid hows at frequencies u ^ k^. As a result, we will need 
to pursue the expansion to higher order in the derivatives of the scalar held than for 
the other helds. Formally we implement this by taking the scalings 

fc ~ e dfj_q ~ dfj,b ~ ~ (2.8) 

and then constructing the solutions as a perturbation series^ in e. 

It is worth emphasising that, beyond leading order, the equilibrium conhguration is 
no longer the translationally invariant black brane (2.2). This is because we can turn 
on position dependent sources for the scalars, = /ca;_ 4 , without inducing a huid 
how. The equilibrium solutions in the presence of these sources can be constructed 
order by order within our derivative expansion. To do this we simply need to set 

= (1, 0, 0, 0) and = d^Jj = 0, but retain the terms that arise from derivatives of 
the scalar helds. The resulting solutions are equivalent to expanding the black branes 
of [11] perturbatively in k. 

Secondly, whilst the huid-gravity correspondence provides us with the full non-linear 
hydrodynamics, we are ultimately interested in linearising these results around equi¬ 
librium. In order to streamline our calculation, we will therefore systematically ignore 
certain terms that will not appear in the linear theory. In particular, we will neglect 
terms of the form 

u^u"' ducl)^^ (2.9) 

^This expansion can be controlled by fixing the ratio /r/T and then taking the large T limit. Working 
perturbatively in derivatives is then valid provided that all fields vary on length scales I ^ T~^. 
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For our choice of scalar sources these terms would be quadratic in the fluid velocity 
and so would play no role in our discussion. Note that in contrast, if we had wanted 
to study forced fluid flows in which was non-zero, then these terms would be 

important and would need to be determined. 

With this caveat in mind, our goal is to pursue this derivative expansion to evaluate 
the constitutive relations for and to 0{e'^) and to calculate (O^) at 0{e^). This 
is enough information to obtain the Ward identities (1.2) at 0{e^) and hence study 
non-trivial aspects of momentum relaxation in the boundary theory. 

Scalar solution at 0{e) 

The hrst thing we wish to do is to determine the scalar prohle at leading order in 
our expansion. The fluid/gravity correspondence works by solving for the correction 
to the ansatz (2.6) ultra-locally - i.e. at a specihc point. More precisely, we use the 
Lorentz invariance of the zeroth order background to choose coordinates in which we 
have = (1, 0, 0, 0) at the point = 0. Note that this new coordinate scheme, in 
which the fluid is at rest, is different from the laboratory frame dehned by the sources 
(j)^^ = kx_A. To emphasise this we use indices i = 1, 2, 3 to refer to the three directions 
orthogonal to the fluid flow. 

We can then proceed by writing the scalar held in terms of the ansatz (2.6) and a 
correction piece 

0^(r, x^) = 0^^ (2.10) 

where the correction piece 0^^ can be determined by expanding the scalar wave equation 

=0 ( 2 . 11 ) 

up to 0{e) and then solving at the point x^ = 0. This results in the differential equation 

( 2 . 12 ) 

which is straightforward to integrate to hnd the solution 

0W(r) = ^ (2.13) 

Note that we have hxed the two constants of integration by demanding regularity 

at the horizon and normalizability in the UV. From this solution we can extract the 
expectation value of the operator, {Oj ), from the UV asymptotics as 

- 167rGjv(0^) = lim y/^g^''dM(t>A + ^^M(\/^7^^^Af0.A) (2.14) 

r—>-oo Z 
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where 7 ^^^ is the induced metric on the boundary. At leading order this yields 


{Oa) = 




IOttGat 


IGttGjv 




( 0 ) 


(2.15) 


Metric and gauge field at 

Our next task is to determine the leading corrections to the metric and gauge field. 
Once again we pick coordinates so that = (1, 0, 0, 0) at = 0. We can then expand 
the ansatz ( 2 . 6 ) around this point to first order in derivatives to get an expression for 
the metric 


9MNdx^ 


dx^ 


2 dvdr — r'^f{r, go, bo)dv'^ + r'^dxidx^ — 2x°‘da(didx'^dr — 

2 x“da/3ir^(l - /(r, go, bo))dx^dv - ^ dv'^ 


where at this order we can define spatial components of the fluid velocity via = 
(l,/9i(a:“)) and we will denote the charge and energy densities at = 0 as go,fco- 


Likewise we can expand the gauge field to get 


Afjdx’^ = - 


73 


rpZ j rpZ 


(2.16) 


Just as in the scalar case these fields are not solutions to the equations of motion. Our 
goal is to find the correction pieces fields 


_ ( 0 ) , ( 1 ) 
9mn — 9mn ' 9mn 

_ 4(0) , 4(1) 


M 


(2.17) 


satisfy the Einstein-Maxell equations 


Gmn — Q 9 MN + 2 


Fmp + ~^9mnFpqF^^ 




M 

N 


_ 1 
“ 2 

= 0 


dM(t>AdN(t>A — ■^9MNdp(j)Ad (pA 


(2.18) 


up to order 0{e^). Note to obtain the equations at this order, it is necessary to use the 
solution for the scalar up to 0{e), which we determined in the previous section. 













Conceptually, the resulting equations fall into two classes. Firstly, a subset of these 
equations are constraints which only act on the boundary data. These are nothing other 
than equations of motion for the currents (1.2) evaluated at 0{e^) . These tell us that 
at any given order in our expansion, the various derivative terms of the hydrodynamics 
helds are not all independent. Rather charge conservation implies 

dvQ = -qodiPi (2.19) 

whilst the equation for the stress tensor is equivalent to 

3dyb = -AhodiPi + 

d,h = (2.20) 

The remaining Einstein-Maxwell equations are a set of coupled ordinary different 
equations that can be used to determine the corrections and . The calculation 
is a straightforward generalisation of the analogous computation in the translationally 
invariant case [23, 24], but with additional terms arising from the scalar stress tensor. 
The details of the equations and expressions for the resulting metric can be found in 
Appendix A. 


Having determined the solution for the metric and the gauge held at leading order, 
it is straightforward to extract the stress tensor in the boundary from the dehnition 




IOttGat 


lim r® 


2{K 




7^"iF) + 67^" +-(V^c/.^V" 




( 2 . 21 ) 


where 7 is the boundary metric, V the associated covariant derivative and is the 
extrinsic curvature. Simlarly the current is given by 

{/<•) = lim (2.22) 

47r(jrAr r-^oo 

Note that due to the ambiguities in dehning the huid velocity and hydrodynamical 
variables out of equilibrium, it is necessary to make a choice of frame in order to obtain 
a unique expression for these currents. In particular we specify our frame by demanding 
that charge and energy densities are held constant and given by bo,qo- We can then 
dehne the huid velocity by demanding that the leading correction to the stress tensor 
obeys 

= 0 (2.23) 

which is commonly referred to as the Landau frame condition. With these specihcations 
we can then read oh the constitutive relations. 
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Thermodynamics 

Although our primary motivation for constructing this derivative expansion was to 
understand the transport properties of these theories, we need to momentarily pause 
and consider their thermodynamics. Once we turn on sources for the scalar fields, the 
thermodynamics of our theory receives corrections. To determine these corrections, we 
set the derivatives of the hydrodynamics fields to zero d^q = d^b = 0 and set the fluid 
to rest = (1, 0, 0, 0) everywhere. 


In the construction of our bulk solutions, we have been working in an ensemble 
with fixed charge density and energy density. The expressions for these quantities are 
therefore unchanged 


3b _ V3q 
IGtiGm ^ AtiGm 


(2.24) 


Conversely the conjugate thermodynamics variables s, T, p and P receive corrections 
that we need to determine. The origin of these changes is that in the presence of sources 
for the scalar fields the (r) component of the metric changes. In particular we have 
that 

gvv{r) = -r'^fir) + = -r^f{r) + (2.25) 


This extra piece shifts the position of the horizon (defined by gwi.'fh) = 0) to 


rh = rQ + 


4ro(4-2g2) +••• 


(2.26) 


where the dots indicate higher order terms in the derivative expansion and Q = q/r^. 
We can see that there are corrections to the entropy density 


rl 1/3 SPro \ 
4^“4^V° + 4(4-2Q2)j + 


(2.27) 


as well as to the chemical potential 


\/3q \/3q y/Sqk'^ 

^ ^ ^ “ 4r^(4 - 2 g 2 ) + • • • 

Similarly we find the temperature is now given by 

4xr = -gUr,) = r„(4 - 2Q=) + + ■ ■ ■ 


(2.28) 


(2.29) 
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and finally pressure is 


P - l,p + sT-e- + + ... (2.30) 

These corrections to the thermodynamics will later be important when we discuss the 
boundary hydrodynamics beyond leading order in k. 


Constitutive relations 


Having extracted the thermodynamics we can now reinstate the other fluctuations 
dfj,q, dfj_b and consider a general fluid velocity u^. At zeroth order in the derivative 
expansion, the stress tensor is simply that of a perfect fluid 




3b 


IGtiGn 




IGttGjv 


piip 


(2.31) 


where + u^u'^ is the projector onto directions perpendicular to the fluid 

velocity. At the correction to the stress tensor (in Landau frame) is given by 




2rl 


167rG 


-a 




N 


32nG 




N 


(2.32) 


where the sources appearing above are 

= P^“P"^a(aM^) - ^P^’^daU^ 
o 

5^305^ (2.33) 

o 

Similarly we can extract the electric current. At leading order this is simply propor¬ 
tional to the fluid velocity but receives subleading corrections proportional to deriva¬ 
tives 




AttGn 


V3 M + 1 
AttGn ^Mtq 


+ 3qu^dxu,) 


1 3\/3 q 

AttGn 3rl M 




(2.34) 


where M = b/vQ. Finally it will be convenient to note that the constraint equations 
(2.20) allow us to eliminate P^'^u^d\Uu in favour of P^'^d^b and We can 

therefore write the constitutive relation for the current in the form 






\/3 M + 1 
dvrGAf 4Mro 


P^’^{d,q 



1 3^3 q^ 

At^Gn 16 




(2.35) 
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Equation (2.35) is our main result so far. The first term two terms, proportional to 
the fluid velocity and derivatives of q and b, are just the usual constitutive relations 
of relativistic hydrodynamics. The novel piece is the additional term proportional to 
the derivatives of the scalar helds. The existence of these subleading corrections in the 
constitutive relation was not accounted for in [2] and will be necessary to describe the 
transport properties of these holographic theories correctly. 


Scalar solution at 0{e^) 

Our hnal task is to calculate the expectation value of the scalars at next order in this 
derivative expansion. The calculation proceeds in entirely the same manner as before. 
Firstly recall that the solution for the scalar at hrst order can be written as 


+ 0 ? 

where we determined the leading correction 

Q, b) = Fi{r, q, b)u>^df,(pA 

with the radial dependence given by the integral 



(2.36) 


(2.37) 


(2.38) 


We now £x coordinates such that = (1, 0, 0, 0) at = 0 and expand the solution to 
0{e^). Note that expanding Fi{r,q,b) introduces a dependence on the derivatives of 
the charge and energy densities. We then search for a correction piece (t)\{r) such that 

4>A = 0^^ + 0^4^ + 0^4^ (2.39) 


solves the scalar wave equation up to terms of O(e^). Note that this involves evaluating 
the wave-equation in the C?(£^) geometry we determined in the last section. 

The resulting equation for 0^^(r) takes the form 

- ('<"))' = ^A{r, diq, d^^i) (2.40) 


where an expression for the source term Sa can be found in Appendix A. Integrating 

this equation twice and imposing the appropriate boundary conditions allows us to 

( 2 ) 

solve for the scalar helds (j)\ . From these solutions we can extract the expectation 
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values using (2.14) as 
{Oa) 



1 QqP'^'^ducj)^^^ {d^q + Squ^dxu^) 

WttGn 2AM Tq 

where iS^ are related to various derivatives of the fluid velocity 


(2.41) 


= {u^dxu>^ ± 


(2.42) 


and A is a complicated function of the dimensionless variable Q. For the neutral theory, 
Q = 0, we have that A = ln2/2 but in general we have only been able to calculate it 
perturbatively (see Appendix A). Fortunately for us, the precise value of A will not be 
important in understanding the physics of the boundary theory. 

In principle we could continue performing this £ expansion to obtain the hydrody¬ 
namics of the boundary theory up to any order in the derivative expansion. However, 
we will not pursue this expansion any further here. Rather, we now have enough infor¬ 
mation to study momentum relaxation, which we address in detail in the next section. 

3 Linear response of the boundary theory 

In the last section we explained how to use the technical machinery of the fluid/gravity 
correspondence to derive the boundary hydrodynamics of our holographic theory in a 
derivative expansion. We now wish to address a particular physical application of these 
results, which is to the transport of charge and heat in the bonndary QFT. 

There has recently been a large amonnt of interest in calculating the thermoelectric 
transport coefficients of holographic models with broken translational invariance. These 
coefficients are dehned in terms of the linear response of the system to electric helds 
and thermal gradients as 



where is the heat cnrrent. Then a is just the usual electrical conduc¬ 

tivity, a is the thermoelectric conductivity and R the heat condnctivity. 
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For many holographic models, exact expressions are known for the DC (cu = 0) 
transport coefficients in terms of horizon data. In particular, for the 5-dimensional 
Einstein-Maxwell-Scalar model (2.1) with translational invariance broken by the sources 
0^ = kxj, these read [11, 15] 


(^DC 

C^DC 

knc 


rp 47rp^ 

AitGn k‘^s 

iirp 

AtisT 


(3.1) 


Importantly these results are exact - that is they are valid for any value of k, includ¬ 
ing when momentum relaxation is very strong. At leading order, that is the 

holographic formulae agree with the hydrodynamic model studied in [2]. In contrast 
at 0{e^) there are differences and the results of [2] no longer apply. Very recently, 
a similar story was found to hold for the low-frequency AC conductivity of the 3-1-1 
dimensional analogue of this theory [19]. 


Given the discussion so far in this paper, such discrepancies should not come as 
a surprise. Rather, we have already emphasised that beyond leading order in the e 
expansion our expressions for the constitutive relations for and (O^) receive 

corrections from the scalar helds that were not present in [2]. Our goal now is to use 
to calculate the thermoelectric conductivities up to 0{e^) using the hydro dynamical 
model we derived in Section 2. We will hnd that, upon taking into account these 
corrections, our results precisely agree with the exact holographic formulae. 


Linearised Hydrodynamics 

The equilibrium conhguration of our holographic model is dehned by placing the 
fluid at rest, = (1,0, 0,0), and setting the gradients d^q, and (9^5 to vanish. The 
sources however, remain non-zero and break the translational invariance of the 
boundary theory. In Section 2 we computed the relevant thermodynamics of this model 
up to 0{k‘^). In order to calculate the thermoelectric transport coefficients we need to 
consider perturbations away from equilibrium 

/i(x, t) ^ p + t) 

T{x,t) ^ T + ST{x,t) (3.2) 

In this paper we will be only interested in studying the frequency dependence of the 
transport coefficients and so consider the perturbations and d^T to be position (i.e. 
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x) independent. We then take a flnid velocity of the form 


u‘‘ = (l,'Uj,((),0,0) 


(3.3) 


We now proceed to linearise the hydrodynamical model derived in Section 2 in these 
pertnrbations and ST{x,t). 

After doing this, the constitntive relation for the momentnm density can be 
written as 


rj-iOx 


4^0 

IGttGn ^ 


(3.4) 


where we have nsed (2.30) and the ... indicate terms of higher order in the e expansion. 
The corresponding expression for the electrical current is 


Jt. 


\/3g 


dvrG 


“'Ut. 


N 


\/3 M + 1 _1 

AttGn 4Mro Ab ^ AttGn 166^ 


KUx + 


pUx (^Q I dxP p 


dxT 


3y/3r 


AttGn 166 ^ 


k^Ux + 


(3.5) 


with (Tg being the ‘universal conductance’ of relativistic hydrodynamics 

_ ro f sT y 

47rGjv[e + Fj 


(3.6) 


To obtain the hnal line in (3.5) we have traded the perturbations dq, 6b with those of 
6fi,6T using 


6q 

6b 



6p + 

T 

6p + 

T 



(3.7) 


We can now explicitly see that the leading (in k) terms in these expressions are sim¬ 
ply the usual constitutive relations for the currents within relativistic hydrodynamics. 
However, the presence of the scalars has introduced extra corrections proportional to 
the fluid velocity which will enter the transport properties at 0{e^). 
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Before proceeding to study momentum relaxation through the Ward identity, it turns 
out to be extremely convenient to rescale the fluid-velocity in order to simplify these 
expressions. We therefore dehne a new fluid velocity via 

Ux = (1 H—+ ■ ■ ■ (3-8) 

so that even up to we can continue to write the momentum density as 

= (e + P)vx ... (3.9) 

The constitutive relation for the current can then be expressed in terms of Vx as 

7 fa , 70 

Jx = pvx - I OxP - 1 + Vx + ... (3.10) 

Finally we need an expression for the heat current Qx = — fiJx. This reads 

^ fa.. 

Qx ^Tvx H“ fi(jQ {oxf^ ^ T / AttT ^^ * * * (3.11) 

where all the thermodynamic factors appearing in these expressions are accurate up 
to 0{e‘^). Written this way the constitutive relations take a remarkably simple form - 
the only novel terms are the 0{k‘^) pieces appearing in the electrical and heat currents 
which can naturally be expressed in terms of thermodynamic quantities. 


Momentum Relaxation 

The constitutive relations for the currents are not enough to dehne our model, we 
also need to supplement them with the equations of motion for the currents 


= 0 

(3.12) 


The second of these equations is simply the Ward identity for momentum non-conservation 
in the presence of the position dependent sources for the scalars [11]. The relax¬ 
ation of our huid velocity Vx to equilibrium is then described by the x component of 
the stress-tensor equation of motion. We can evaluate this for our perturbations using 
the expression (2.41) for (Oa). At order 0{e^) this equation reads 


+ dxP 




dvr 


■Vx + ... 


k^s 


j 


47r(e -I- P) 


(3.13) 
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which is precisely the theory studied in [2] where the momentum relaxation rate is 
identified as 


T 


-1 


k'^s 

4ir(e + P) ^ ■ 


(3.14) 


However, just as in the constitutive relations for the currents, there are subleading 
corrections to (3.13) that differ from [2]. Although our expression for the scalar ex¬ 
pectation value (2.41) is rather complicated, it simplifies dramatically for linearised 
perturbations. Evaluating the Ward identity gives at 0{e^) gives the remarkably suc¬ 
cinct equation‘s 


dtT^^ + = 



1 ( 

SVx+ rj. 1 

An 

k'^ Qx Xk'^VQ 

An 

T IGttGat 




/i 


T 

dtVx + ... 


dxT 


47rT2 


Xk'^r'^ 

'0 o I 

1 C n + ... 

IdttG jv 


(3.15) 


where we have noted the quantity in square brackets is nothing other than the entropy 
current Qx/T. 


Once we go beyond leading order in the e-expansion we therefore find that the right 
hand side of (3.15) is not simply proportional to the momentum density Surpris¬ 
ingly it is the heat current, which differs from the fluid velocity through derivatives 
of the hydrodynamical fields, that more naturally appears in the equation for momen¬ 
tum relaxation. It is not clear to us whether there is any deep reason why this should 
be the case, but we will shortly see that it is closely related to the existence of simple 
formulae (3.1) for the DC transport coefficients of this model. 


Thermoelectric response coefficients 

We can now turn to our final question of interest which is to extract the response 
coefficients from our linearised hydrodynamics. To do this we use the Ward identity 
(3.15) to get an expression for the fluid velocity in terms of the perturbations dxH and 
dxT. Using the thermodynamic relationship 5P = s6T + p6fi we find 

1 /dvru ucTo U^Pcro\ ... 1 dvr, 


Vx=- 


lUJT 




k'^s 


sT 


g2T2 


lUT 


k^ 


dxT + ... (3.16) 


where the momentum relaxation rate, r ^ is determined up to 0{e'^) to be 


r ^ = 


sk'^ 


1 - 


p^aqk'^ 


\k‘^rl 


47r(e + P) [ -|- P) 


+ 


(3.17) 


"‘In deriving this equation we found it helpful to recall (2.30) to write 5Txx = A:^roi5ro/(167rG7v) 

and again make use of the constraints (2.20). 
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Since the corrections to r ^ depend on A they are some non-trivial fnnction of Q that 
we cannot in general express analytically®. 


By inserting (3.17) into the constitntive relations for the cnrrents we can calculate 
the thermoelectric transport coefficients perturbatively® in ca ~ ~ At we 

have the result 


a{u)) 

R,(u) 


1 

1 — iuT k?s 
1 dvrp 

1 — iujT k'^ 

1 dvrsT 

1 — iujT 


(3.18) 


At this order one simply sees a ‘coherent’ metal, i.e. we have a well-dehned Drude 
peak at low frequencies. As expected these expressions match the results of [2] with 
the momentum relaxation rate (3.14). 


However the corrections to the constitutive relations and the Ward identity enter at 
0{e^) and the transport coefficients are no longer equivalent to those in [2]. Rather we 
hnd 


a{u) 

a{u) 

R,(u) 


(^Q + 


ag + 

Kq + 


1 

1 — iUT 
1 

1 — ioJT 
1 

1 — iUT 


4:71 

k‘^s 

dvrp 

/ 47rsT 


+ (To — (Jq 


«(? 




(3.19) 


where we have dehned ag = —p/Tag and Rg = p^/Tag and (Tq = ro/47rG)v. In the 
DC limit cj ^ 0 these expressions reduce to the exact formulae in terms of horizon 
data (3.1). 

The AC conductivity can be seen to be composed of two distinct pieces. In addition to 
the coherent Drude peak, there is now an extra ‘incoherent’ (i.e. frequency independent) 
contribution. This incoherent contribution to the transport coefficients simply arises 

^Note that a similar expression for the momentum relaxation rate in the 3+ld analogue of our 
model was presented in [19]. We have checked that the lower dimensional equivalent of (3.17) agrees 
with their expression. 

®For our perturbations, this is simply the usual Kadanoff-Martin prescription for extracting hydro- 
dynamic correlation functions [2, 27]. Note, however, that for spatially modulated perturbations then 
Sfi{x),ST(x) have their own dynamics and the computation of transport coefficients is more involved. 
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from the parts of the constitutive relations that are proportional to the derivatives 
dxT. Conversely, the coherent piece corresponds to the part of the currents that 
is proportional to the fluid velocity and so relaxes at the rate 

Reassuringly, an identical structure to (3.19) was recently found in the 3+1 dimen¬ 
sional analogue of (2.1) through an explicit computation of the Green’s functions of the 
boundary theory [19]. The fact that we precisely agree with these results is a strong 
check on the consistency of our hydrodynamical description. 


One of the most intriguing consequences of (3.19) is that it reproduces the DC 
conductivity in a rather non-trivial way. In particular, the 0{e^) pieces in (3.1) arise as 
a result of rather surprising cancellations between the incoherent and coherent pieces. 
It was not clear from the discussion in [19] why the sum of these pieces should be 
simple. Our hydrodynamical formulation allows us to see that these cancellations are 
closely related to the form of momentum relaxation in these models. To see this we 
should look more carefully at the Ward identity 


dtT^^ + d^P 


^0. 

47r T 


Xk'^rl 

IQtiGm 


dtVx + ■ • • 


(3.20) 


In the DC limit, we can ignore the time derivatives in this expression and simply balance 
the forces. Consequently we must satisfy 


Qx 


AttT 


dxP 



( dvrsT 


dxT 


(3.21) 


from which we can immediately read off the expressions for a^c and Rdc' ■ The sim¬ 
plicity of these DC transport coefficients therefore arises from the fact that it is the heat 
current, not the momentum, which appears on the right hand side (3.20). That the 
hnite uj results are more complicated is then not surprising, since the time derivatives 
in (3.20) now mix the momentum and heat currents together. As a hnal observation 
it is natural to speculate that, since the holographic formulae (3.1) are exact, an equa¬ 
tion like (3.20) involving the heat current should continue to hold at all orders in the 
derivative expansion. 


4 Discussion 

In this paper we have used the fluid/gravity correspondence to construct a hydro¬ 
dynamical description of a simple holographic theory in which translational invariance 

more detailed study of the DC limit in these models was subsequently performed in [28]. 
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is broken by scalar fields. We systematically derived the constitutive relations for the 
boundary current, stress tensor and expectation value of these scalars in a derivative 
expansion. Beyond leading order in this expansion, we found that the resulting hydro¬ 
dynamics was different to that studied in [2]. 

We then looked at charge and heat transport in the boundary theory by linearising 
this hydrodynamics about equilibrium. This allowed us to calculate the low-frequency 
thermoelectric response coefficients perturbatively in the strength of momentum re¬ 
laxation. Our model reproduced the known formulae for the DC conductivities in 
terms of horizon data. We also found that the AC conductivity could be decomposed 
into ‘coherent’ and ‘incoherent’ parts in the same manner as was recently seen in [19]. 
In summary we have constructed, for the first time, a hydro dynamical description of 
holographic charge transport with momentum relaxation. 

Now that we have developed this model, there is a myriad of possible directions 
for future research. In this paper we have only considered the frequency dependence 
of the linear response functions. There is much more information contained in the 
spatially resolved transport coefficients. It would be interesting to develop a detailed 
understanding of diffusive processes in these theories, analogous to the discussion in 
[26]. Secondly, since our description differs from that presented in [2], it is natural 
to revisit the motivating questions of that work. Through applying the fiuid/gravity 
correspondence in a magnetic field we could construct a theory of magnetotransport 
[29-32] and study the frequency dependence of the Hall angle [33]. 

Finally, in this paper we have focused on the simplest holographic model in which 
we could address the question of momentum relaxation. Nevertheless the techniques 
developed in Section 2 can be extended to more general settings. In particular, although 
we expect broadly similar results, the correspondence should be developed for 3-1-1 
dimensional bulk actions [34]. It would also be particularly interesting to consider 
more general, in particular inhomogeneous [4, 5, 10, 35], sources for the scalars and 
understand how the resulting hydrodynamics differs from that presented here. 

We hope that developing such hydrodynamical descriptions will provide a physical 
picture of transport in holography and hence allow for a greater understanding of the 
potential implications of these models for real-world systems. 
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A Details of the fluid/gravity calculations 

In this appendix we provide, for the beneht of the interested reader, some further 
details about the fluid/gravity computations we used to obtain the constitutive relations 
presented in Section 2. We have already explained, in the main text, how to compute 
the scalar solution to 0{e) and so begin with the leading corrections to the metric and 
gauge held. 

Metric and gauge field at first order 

Our goal in this section is to solve the Einstein-Maxwell equations (2.18) at = 0 
to determine the 0{e‘^) corrections to the the metric and gauge held As is 

standard in these calculations, we proceed by decomposing these equations into their 
scalar, vector and tensor pieces. 


Scalar modes 

We begin, following the discussion in [23], by parameterising the scalar parts of the 
metric and gauge held as 


where gu and g^r 


9ii^ = r^hi(r) (no sum) 
= ki{r)/r‘^ 

9vr = -^^i(^) 

(1) ^ VSwijr) 
y 2r2 

are related by a choice of gauge 


(A.l) 


Tr 




0 grr = 0 gr/j. oc = 0 


(A.2) 
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As we mentioned in the main text, the Einstein-Maxwell equations can be decom¬ 
posed into constraint equations and dynamical equations. The constraint equations 
are dehned as those that arise from the dot product of the Einstein-Maxwell equations 
with the vector dual of dr. 

Importantly, a subset of these constraint equations correspond to the hydrodynamical 
equations of motion of the boundary theory. The hrst such constraint equation, arising 
from the Maxwell equations Mjv = 0 is given by 

g^^Mr + = 0 (A.3) 

which evaluates to the requirement that we have 

dvQ = -qodijdi (A.4) 

and is equivalent to the equation for conservation of the boundary current 

= 0 (A.5) 

Similarly we have the constraint equations arising from the Einstein equations Emn = 
0. Firstly we can take 


+ = 0 (A.6) 

which reads 

?>d^h = -4hodi(3i + (A.7) 

and corresponds to the scalar component of the Ward-identity for momentum relaxation 

(A.8) 

Finally there is one more constraint equation arising from 

g^'^E„ + g^'^E,r = Q (A.9) 

This equation does not act on the boundary data, but rather imposes a relationship 
between hi(r), ki{r) and wi{r) 

2dil3ir^ + 12r®hi(r) -|- 4qoWi{r) — bor^h[{r) + 3r'^h[{r) — r^k[{r) — 2qorw[{r) = 

(A.IO) 
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Together with these constraint equations we also have the remaining dynamical Einstein- 
Maxwell equations. In particular the rr-component of the Einstein equation reads 

3h"(r) + lhh\{r)/r = —(A.11) 

whilst the n-component of the Maxwell equation is 

rw'l{r) — w'^{r) — 6 gohi(r) = 0 (A. 12 ) 

These three equations ((A. 10) (A.ll) and (A. 12 )) can now be solved for the undeter¬ 
mined functions hi(r), fci(r), tci(r). 

In general solving these coupled equations results in rather ugly expressions for the 
metric. However, the origin of many of these complications are the terms proportional 
to 4>^^'{rY As we described in Section 2 , these will not contribute to 

the the linearised hydrodynamics with our choice of sources (2.7). In order to make our 
discussion more presentable, we will therefore neglect these terms in equations (A. 10 ) 
and (A.ll). 

If we do this we can then hrst solve for hi{r) using (A.ll) to get 

hi(r) = C'i/r^ + C '2 (A. 13) 

Plugging this into (A. 12) yields wi{r) 

wi{r) = + C 4 - goCi/r"^ (A. 14) 

Finally, we can then use (A.10) to solve for fci(r) 

ki{r) = + C 5 - 2 qo/r‘^C 4 + { 2 ,ql/r^ - hQ/r'^)Ci (A.15) 

All that remains is to hx the various constants of integration that appear in the solution. 
Firstly, the constants C 2 and correspond to non-normalizable modes and so can be 
set to zero in the UV. Secondly, the constants C 4 and C 5 can be absorbed into the 
dehnition of charge and energy density and can also be ignored. Finally the constant 
Cl can be removed by a redehnition of the radial coordinate r r(l -|- A/C). As a 
result, we can conclude that the solution in the scalar sector is 

hi{r) = wi{r) = 0 ki{r) = ‘^CdiC + ^’^^(^*0®)^ (A.16) 

where once again we emphasise that in general there are terms proportional to (9^0®)^ 
in these functions which we are neglecting for the sake of simplicity. 
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Vector modes 


We can now proceed to calculate the corrections to the vector channel, that is we wish 
to solve for the corrections This will be particularly important for us because 

it is this calculation that determines the constitutive relation for the electrical current. 
In this sector we have one constraint equation, arising from 

g^^Er^ + g^^E,i = 0 (A. 17) 


which gives 

dib = -Abod^Pi - rld^cj)^^ (A.18) 

and is equivalent to the vector component of the Ward identity for momentum relax¬ 
ation. The corrections to the metric and gauge held are determined by the dynamical 
equations. In particular the Einstein equations Eri read 

+ 4\/3goaf (A.19) 


whilst the integrated Maxwell equations yield 

A/(r)oy' + 71^ + ^{qAP, + 9(9)) = A’ (A.20) 

where Dj is a constant of integration that we will eventually determine by requiring 
regularity at the horizon. These equations can be solved by using (A. 19) to write 

-t- rV?'(r’)di05^ + 4rgl2 - r^N^P')') (A.21) 

4v 3go 

and inserting this into the Maxwell equation to get a second order differential equation 
for giy 

r^f{r){4:rgiy - r^irgiy')') + Uq^gi^/r^ = - 6qo{diq + 3qod^(3i)/r 

+ iVSqoDl - 

-h {3r^f{r) + 12ql/r)dyl3i (A.22) 


which has the solution 

giy = dv^ir + 


= r’V(r) 


VI 

r*00 


dr 


Ai 

/(r )2 V 3 r® 


2V3qoDl 


+ 


+ 


2P 


a(o) I 


5.07 ^* 0 ) 4 ' + 


+ A" 


where A* = Qqo{diq + 3qod^/3i). 
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Once again we need to fix some constants of integration. Df corresponds to a non- 
normalizable mode and hence shonld be set to zero. Upon calculating the stress tensor 
of the dual theory one hnds that a non-zero Df gives a contribution to the momentum 
oc Df. In Landau frame = 0 this means we must set Df = 0. Finally D] is 

hxed by requiring regularity of at the horizon. From (A.20) this requires that 

+ d,q) -Dl = Q (A.23) 

at the horizon. Using the solution (A.25) we hnd that this constraint implies that we 
take 


D] 


\/3 / Mq + 1 \ 
4ro \ Mq ) 


{dig + 3qody(3i) + 


3\/3 go 
8r^ Mo 




(A.24) 


where Mq = bo/r^ is the dimensionless mass of the black brane. This determines the 
metric correction to be 




Mq + I 
4Moror'^ 


Gqoidiq + SqodyjSi) 


r® 2P AMorff'^ 


di(j) 


( 0 ) 

A 


(A.25) 


and we will leave the gauge held as implicitly determined by this result through (A.20) 
and the requirement of normalizability in the UV. 


Tensor modes 

Finally we can determine the solution in the tensor sector. Here we parameterise the 
tensor components of the metric via Qij = where aij is a tensor of SO{3). In this 

sector we only have a dynamical Einstein equation which reads 


(rV(r)ak)' = -r(j)ij - 6r^aij 

(A.26) 

where we have dehned the source terms 


(405^)2 

(A.27) 

_ dil3j + djl3i ^ ^ 

(A.28) 


Integrating this equation and imposing regularity and normalizability gives the solution 


~ 2 


coo :^2 

df ° 




2a, 




coo 3 

df ° 


y(y)yE 


(A.29) 
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which completes the determination of the metric. The expectation values of the stress 
tensor and current at can then be read off from equations (2.21) and (2.22). Co- 
variantizing these expressions in terms of the fluid dynamical variables q{x),b{x) and 
u^{x) gives the constitutive relations presented in the main text. 


Scalar field at second order 

We can now present the details of how we determine the scalar field solution to 
0{e^). Expanding the scalar wave equation about the point = 0 and using the 
metric determined in the last section yields the following differential equation for the 
correction® </>^^(r) 

_ (r®/(r)0^^^(r))' = (r + 2r®F((r) + 3r^Fi(r))S'^ — 2rS^ 

(A.30) 


where Fi{r) is the radial profile 




POO 


Fi(r) = / dr 


f5/(r) 


(A.31) 


and F 2 {r) is defined by 


rVg(r^ + rro + rg) 


F ^_ ' ' ov' ^ ' 0 ^ ' o; _ 

^ 12(r + ro)(-go + + ^o)) 

The source terms are related to derivatives of the fluid velocity 


= dyf3idi(p^A^ ± ^di(3id^(p^A^ 
g we licive eiiOSeii lu eiiimiiciie aiiy leiiiis p 

to dib using the constraint (A. 18). It is straightforward to integrate (A.30) 
- (rV(r)02^'(r)) = 


(A.32) 


(A.33) 


and we note that in deriving (A.30) we have chosen to eliminate any terms proportional 
constraint (A. 18). It is straightforward to integrate (A.30) to obtain 


^V'(r\\= I df(r + 2r®F((r) + 3r®Fi(r))S'^ — f dr2rS'^ 

J rc\ 


/ro 


•J TO 

- rh^y^di(ty^ - E^a 


+ F2{r 

(A.34) 

®Note that the ignored pieces ~ in ki{r) and hi{r) would contribute additional terms to 

this equation which we are continue to neglect. 
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where the constant of integration is hxed by demanding regularity at the horizon 
to be 



ro 


8(2 - Q2 




3ro 

8A4q 



2AMo4 


(A.35) 


with Qo = qo/r^. The full solution for 4>j^\r) can then be obtained by a further 
integration and requiring normalizability in the UV. Finally we can extract the {OA) 
aX = 0 as 


167iGn{0_a) 


-r^d^A - ^^o^A + ^ 0 *^. 


2o- 


( 0 ) 


-rodv(pA 


ro 


O^A 

- >^rlS+ + 


E 


( 1 ) 

A 


-Qt 

di4)^^\d^4)Pdi(l)P] 


8(2 
3ro 
8Mc 

Qqodi(j)A\9iq + 8qody(5i) 
2AM,rt 


(A.36) 


Covariantizing this expression gives the constitutive relation (2.41). As we mentioned 
in the main text we have not been able to obtain an exact expression for A. It is 
possible, although not very illuminating, to calculate A perturbatively, for which we 
hnd the leading terms 


A 


ln2 


n2 

+ (3-ln4)^ + ... 


(A.37) 


Fortunately, the exact value of A did not play an important role in our discussion of 
transport in Section 3. 
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